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The electronic structure of a monolayer of single-molecule magnets Mni2 on a Au(lll) surface is 
studied using spin-polarized density-functional theory. The Mni2 molecules are oriented such that 
the magnetic easy axis is normal to the surface, and the terminating ligands in the Mni2 are replaced 
by thiol groups (-SH) where the H atoms are lost upon adsorption onto the surface. This sulfur- 
terminated Mni2 molecule has a total magnetic moment of 18 \ib in the ground state, in contrast 
to 20/xs for the standard Mni2. The Mni2 molecular orbitals broaden due to the interaction of 
the molecule with the gold surface and the broadening is of the order of 0.1 eV. It is an order of 
magnitude less than the single-electron charging energy of the molecule so the molecule is weakly 
bonded to the surface. Only electrons with majority spin can be transferred from the surface to 
the sulfur-terminated Mni2 since the gold Fermi level is well above the majority lowest unoccupied 
molecular orbital (LUMO) but below the minority LUMO. The amount of the charge transfer is 
calculated to be 1.23 electrons from a one-dimensional charge density difference between the sulfur- 
terminated Mni2 on the gold surface and the sulfur-terminated Mni2, dominated by the tail in the 
electronic distribution of the gold surface. A calculation of a level shift upon charging provides 
0.28 electrons being transferred. The majority of the charge transfer occurs at the sulfur, carbon, 
and oxygen atoms close to the surface. The total magnetic moment also changes from 18 /j,b to 
20 \ib, which is due to rearrangements of the magnetic moments on the sulfur atoms and Mn atoms 
upon adsorption onto the surface. The magnetic anisotropy barrier is computed including spin- 
orbit interaction self-consistently in density-functional theory. The barrier for the Mni2 on the gold 
surface decreases by 6 K in comparison to that for an isolated Mni2 molecule. 

PACS numbers: 75.50.Xx, 75.70.Ak, 75.30.Gw, 71. 15. Mb 



I. INTRODUCTION 

A nanoscale single-molecule magnet (SMM) comprises 
a few transition metal ions interacting through organic 
and/or inorganic ligands instead of direct exchange in- 
teractions. To reverse the magnetic moment of a SMM, 
a large energy barrier needs to be overcome. Experi- 
ments on bulk forms of SMMs exhibited quantum tunnel- 
ing between different directions of magnetic moments,-^ 
quantum interference between spin paths, 3 and long spin 
dephasing time T^X These properties of SMMs pro- 
pelled interest in utilizing SMMs as information storage 
devices^ spin-based devices^ or materials for quantum 
computation J- A great amount of experiments were car- 
ried out on deposition of SMMs Mni2 or its derivatives 
on gold&^ i 10 ' 11 and silico n 12 ! 13 ! 14 surfaces or on bridging 
them between gold electrode a 15 ! 16 ' 17 . Among thousands 
of synthesized SMMs, [Mni 2 Oi2(CH3COO)i 6 (H 2 0) 4 ] 
(referred to as Mni 2 ) was widely studied due to its large 
magnetic anisotropy barrier (MAB) or magnetization re- 
versal barrier of 65 KJ£ Largely, Mni 2 molecules were 
deposited onto a surface or bridged between electrodes 
in two different manners: (i) through attractive van 
der Waals forces between the surface and Mni 2 with- 
out surface-binding ligands or (ii) via ligand exchange 
with the Mni 2 molecules. In the latter case, for exam- 
ple, carboxylate terminated alkancthiolates, C„H 2 „_|_2S, 
would be surface-binding ligands between the molecules 
and the surface^ or 16 acetate (0 2 C-CH 3 or Ac) ligands 
within Mni 2 could be replaced by ligands 2 C-C6H4-SAc 



that create a direct strong bond to goldJ£ 

Scanning tunneling microscope (STM) image d 9 ' 10 ! 19 
and atomic force microscope (AFM) images 1 ^ on the 
monolayers of Mni 2 molecules revealed that the Mni 2 
molecules in the monolayers are individually distinguish- 
able rather than aggregated on the surface. Photoe- 
mission spectroscopy experiments on the monolayers of 
Mni 2 derivatives showed that the Mn 3d partial density 
of states in valence bands for the monolayers is compa- 
rable to that for bulk Mni 2 i 20 i 21 A low concentration 
of molecules on surfaces makes it challenging to accu- 
rately measure the magnetic properties of a single mono- 
layer of Mni 2 . A recent magnetic measurement on thick 
monolayers of Mni 2 molecules exhibited qualitatively dif- 
ferent magnetic properties from bulk Mn^ji 1 - Very re- 
cently, local magnetic properties of a Mni 2 monolayer 
were measured using depth-controlled /3-detected nuclear 
magnetic resonance (NMR) . This measurement also sug- 
gested that the magnetic properties differed from those 
of bulk Mni 2 »i 4 - In all of the experimental systems dis- 
cussed little was known about the interface between the 
Mni 2 molecules and the surface or electrodes. For ex- 
ample, the following questions have not been answered: 

(i) how the Mni 2 molecules are oriented at the interface, 

(ii) whether the Mni 2 molecules remain chemically in- 
tact at the interface, and (iii) how strongly the molecules 
are coupled to a surface or electrodes. Therefore, it is 
still a controversy whether the electronic and magnetic 
properties of SMMs change due to the interaction with a 
surface. 
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FIG. 1: (Color online) Schematic (a) top and (b) side views of 
a Mni2 molecule adsorbed onto a gold slab via thiol groups. 
The three different Mn ions (regarding symmetry) are denoted 
as I, II, and III in (a). The side view of the six gold monolayers 
with 36 surface atoms per monolayer are shown in (b). In 
supercell calculations, a vacuum layer of 10 A is placed above 
the Mni2 molecule. 



To investigate mainly electronic transport and spin 
filtering through SMMs, theoretical models 2 ^ 2 ^! 2 ^ 
were proposed based on the assumption that individ- 
ual SMMs were not strongly coupled to electrodes, in 
other words, the molecular orbitals remain sharp despite 
their interactions with metal electrodes. Then the An- 
derson Hamiltonian and an effective spin Hamiltonian 
were used with a priori microscopic parameter values that 
need to be determined from atomic-scale simulations. To 
the best of our knowledge, SMMs deposited on a surface 
or bridged between electrodes, have not been yet stud- 
ied using first-principles methods. We investigate, within 
density- functional theory (DFT), the effect of the inter- 
face on the electronic and magnetic properties of Mni2 
molecules on a gold surface, given a particular orientation 
of Mni2 molecules relative to the surface. More specifi- 
cally we obtain the strength of coupling of the molecules 
to the surface and other parameter values that could be 
used in a model Hamiltonian. 

We consider a monolayer of Mni2 molecules adsorbed 
onto a Au(lll) surface through thiol groups as shown in 
Fig. [TJ Here we use the shortest chemical link between 
the Mni2 molecules and the gold surface. Although this 
link may seem too short to create stable Mni2 mono- 



layers in experiments, it has advantages for theoretical 
studies. The link would considerably reduce computa- 
tional cost by requiring a smaller unit cell without alter- 
ing the physical and chemical properties of the system. 
If the molecules can be strongly coupled to the surface, 
this link would be a good candidate. Our DFT calcula- 
tions show that the Mni2 molecular orbitals moderately 
broaden for the short link due to the interaction between 
the molecules and the surface. Electronic charge is trans- 
ferred from the surface to the molecules and the total 
magnetic moment is also modified upon adsorption onto 
the surface. Our model and methods for studying a Mni2 
monolayer on the gold surface are presented in Sec. II. 
The electronic structure and magnetic properties on the 
system are discussed in Sec. III. The conclusion follows. 



II. MODEL AND METHODS 

Spin-polarized DFT calculations are performed with 
plane waves as basis sets, within the Perdew-Burke- 
Ernzerhof (PBE) generalized-gradient approximation 
(GGA), 27 using the Vienna Ab-initio Simulation Pack- 
age (VASP) i 28 i 29 Projector-augmented- wave (PAW) 
pseudopotential a 30 i 31 are employed to take into account 
spin-orbit interaction (SOI) that induces large magnetic 
anisotropy in a Mni2 molecule. Valence electrons con- 
sidered in the PAW pseudopotentials for each atom are 
shown in Table [TJ 5d 10 6s 1 orbitals are treated as valence 
states for Au, and 3p 6 4s 2 3d 5 orbitals are used for Mn. 
Hard pseudopotentials are used for C and O. All DFT cal- 
culations in this study are performed self-consistently us- 
ing VASP (unless stated otherwise) until the total energy 
converges down to l.Ox 1CP 5 eV. To simulate a monolayer 
of Mni2 molecules on a Au(lll) surface (Fig. [JJ), we con- 
sider a monoclinic unit cell of 17.71x17.71x34.00 A 3 in 
which a simplified form of a M1112 molecule (explained 
specifically in Sec. II. B) is attached to a gold slab with 
six monolayers via thiol groups, and a vacuum layer of 
10 A is added above the Mni2 molecule. In this geometry 
the magnetic easy axis of the M1112 molecules (z axis) is 
oriented normal to the surface. We first optimize a gold 
slab and a Mni2 molecule geometry, separately, and then 
combine them to create what we call the 'whole structure' 
illustrated in Fig. [JJ 



A. Gold slab 

A DFT calculation is performed on bulk gold and con- 
vergence of the total energy is checked as a function of the 
number of fc-points, fast Fourier transform (FFT) mesh, 
an energy cutoff for plane waves, and a cutoff in augmen- 
tation charges. With an energy cutoff of 260 eV and a 
cutoff in augmentation charges of 357 eV, the calculated 
equilibrium lattice constant is 4.175 A that differs from 
the experimental valued by 2.8%. In constructing the 
gold slab, in-plane separations among gold atoms are set 
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FIG. 2: (Color online) (a) Geometry of the sulfur-terminated 
Mni2, geometry 2, and its total energy as a function of to- 
tal magnetic moment, (b) Geometry of the thiol-terminated 
Mni2, geometry 3, and its total energy as a function of total 
magnetic moment. 

to be 2.952 A, obtained from the equilibrium geometry of 
bulk gold. With the in-plane separations fixed, the dis- 
tance between monolayers is determined by atomic force 
relaxation until the magnitude of all force components is 
less than 0.01 eV/A. To fully cover a Mni2 molecule on 
a Au(lll) surface, a 6x6x1 real-space supercell (36 sur- 
face Au atoms per monolayer) is used. The surface area 
covered by the molecule is stressed by the dotted circle 
in Fig. [T]Ja). In addition, a gold slab should be thick 
enough to completely screen the Mni2 molecule includ- 
ing thiol groups with height of 12.12 A. This condition 
is satisfied with six gold monolayers 12.06 A high. In 
the whole structure we include the optimized geometry 
for a six-monolayer gold slab with 36 surface atoms per 
monolayer that contains a total of 216 gold atoms. 



B. Isolated Mni2: geometries 1 and 2 

To reduce computational cost, a simplified form of a 
Mni2 molecule, [Mni20i2(HCOO)ig] (denoted geometry 
1), is used in our study. The structure of geometry 1 is 
the same as what is shown in FigJ^a) except that the S 
atoms are replaced by H. The simplified Mni2 differs from 
the synthesized Mni2 molecule in that the 16 acetate lig- 
ands were replaced by 16 formates (HCOO) and the four 
water molecules were removed because they may be lost 



TABLE I: Valence electrons, energy cutoffs, and atomic 
sphere radii for the atoms where the PAW pseudopotentials 
are employed. The given atomic radii are used to calculate 
atomically resolved magnetic moments shown in Table HT1 

Atomic Valence Electron Suggested Energy Atomic Sphere 



Species 


Configuration 


Cutoff Range (eV) 


Radii (A) 


Au 




170-230 


1.50 


Mn 


3p 6 4s 2 3rf 5 


202-270 


1.32 





2s 2 2p 4 


500-700 


0.74 


c 


2s 2 2p 2 


500-700 


0.86 


s 


3s 2 Sp 4 


302-402 


0.95 


H 


Is 1 


200-250 


0.37 



via ligand-exchange during adsorption onto a surface, 8 
Geometry 1 is placed in a unit cell of 25.0x25.0x20.0 A 3 , 
where the magnetic easy axis is along the z axis. Then 
geometry 1 is relaxed with a fixed total magnetic moment 
of 20 hb (Bohr magneton) and an energy cutoff of 600 eV 
until the magnitude of all force components becomes less 
than 0.08 eV/A. Using S4 symmetry of geometry 1, the 
twelve Mn ions are categorized into three classes marked 
as I, II, and III in Fig. QJa). The four inner Mn ions 
(Mn 4+ ) belong to symmetry class I, while the two sets of 
the four Mn ions (Mn 3+ ) in the outer ring belong to sym- 
metry classes II and III, respectively. The spin moments 
of the Mn 4+ ions are antiparallel to those of the Mn 3+ 
ions in the ground state. Our noncollinear calculation 
shows that a collinear spin approximation is good for ge- 
ometry 1. The simplification made for the Mni2 molecule 
does not significantly change the electronic and magnetic 
properties of the molecule. The gap between the major- 
ity highest occupied molecular orbital (HOMO) and the 
majority lowest unoccupied molecular orbital (LUMO) 
is calculated to be 0.24 eV, while the minority HOMO- 
LUMO gap is 1.95 eV. 

To attach the Mni2 molecules onto a gold surface, the 
H atoms in the two formates in geometry 1 are exchanged 
with surface-binding ligands, SH. Then the H atoms in 
SH are removed for the S atoms to be directly bonded 
to the gold surface. This determines the geometry 2, as 
shown in Fig. Ela). The distance between the S atom 
and each bonding C atom is 1.88 A, which was obtained 
from geometry optimization of Mni2-alkane-S-Aui3. The 
Fermi level of the gold slab as well as molecular orbitals 
in its vicinity including HOMO and LUMO are depicted 
in Fig. [3] for geometries 1 and 2 and the whole struc- 
ture. In geometry 2 the S atoms are highly reactive and 
their presence reduces both the majority and minority 
HOMO-LUMO gaps compared to geometry 1. Due to 
the lost H atoms, the HOMO and LUMO of geometry 2 
are shifted downward compared to geometry 1 and the 
total magnetic moment of geometry 2 decreases to 18 /is. 
The calculated total energy of geometry 2 as a function 
of total magnetic moment is shown in Fig. [Ha). Mag- 
netic moments for individual atoms are calculated over 
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TABLE II: Atomically resolved magnetic moments for geometry 2 and atomically resolved difference of the magnetic moments 
between the whole structure and geometry 2 in units of Bohr magneton. The positive sign in the magnetic moment change 
denotes an increase in the magnetic moment upon adsorption. The three different Mn ions in terms of symmetry, Mn(I), 
Mn(II), and Mn(III), are depicted in Fig. 0Ja). Significant increases in the magnetic moments occur at the atoms whose 
magnetic moments are bold-faced. 



Atomic species 


Initial magnetic moment 


Change of magnetic moment upon adsorption 


Mn (I) 
Mn (II) 
Mn (III) 


-2.621,-2.620,-2.573,-2.571 
3.453,3.461,3.530,3.516 
3.541,3.543,3.527,3.534 


+0.019, +0.021, +0.002, +0.001 
+O.O63,+O.O6O,+0.004,+0.004 

+0.002,+0.003,+0.010,+0.010 


O closest to Mn and S 
O (total) 


-0.011,-0.011,-0.028,-0.027 
0.328 


+0.015, +0.014, +0.036, +0.035 

+0.105 


C closest to S 
C (total) 


0.007,0.007 
0.138 


-0.008,-0.008 
-0.016 


S 


-0.236, -0.235 


+0.237,+0.238 



-4.5 
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FIG. 3: (Color online) Molecular orbitals in the vicinity of 
the gold Fermi level (dashed line) for geometries 1 and 2 
and the whole structure. The orbital energies are written for 
the majority and minority HOMO and LUMO. The orbitals 
shown for the whole structure are projected ones onto all the 
Mn atoms only. 



atom-centered spheres with radii given in Table I. See Ta- 
ble II. The S atoms bear substantial magnetic moments 
of 0.24 [is parallel to the inner Mn(I) ions. The magnetic 
moments of the Mn(I) and Mn(II) ions (denoted bold in 
Table II) connected to the S atoms via the O and C atoms 
are considerably smaller than those for geometry 1. 



C. Whole structure: geometry 2 on Au(lll) 

Top and side views of the whole structure used in our 
DFT calculation are shown in Fig. [TJ Each sulfur atom 
is bonded to the closest hollow (face-centered cubic) site 
on the gold slab which is the most favorable configura- 
tion. Attaching a sulfur atom at a bridge site requires a 
0.30 eV higher energy than at a hollow site. This result 
agrees with a calculation on alkane thiols on a gold sur- 
face, for which a hollow site was favored by 0.40 eVi22 
A top site on a gold slab is the most expensive energeti- 
cally. The distance from the sulfur atoms to the gold slab 
is optimized such that both sulfur atoms are as close as 
possible to hollow sites. The distances between the sul- 
fur and the gold atoms range between 2.52 and 2.74 A, 
where the longest bond is due to a slight mismatch of 
placing one of the sulfur atoms at the closest hollow site 
on the gold slab. 

In the whole structure the closest distance between hy- 
drogen atoms on periodic images of the Mni2 (geometry 
2) is 3.35 A. The large intermolecular separation pre- 
vents different Mni2 molecules from interacting with one 
another except for dipolar interactions. In the present 
study the dipolar interactions are not considered. The 
whole structure has a total of 304 atoms and 2886 va- 
lence electrons. Due to the presence of Mn atoms, an en- 
ergy cutoff of 600 eV for plane waves and a cutoff in the 
augmentation charges of 800 eV are used. For the whole 
structure, we sample 2x2x1 fc-points, including the T 
point. Notice that we used a 6x6x1 real-space supercell 
per gold monolayer. Thus, on the plane parallel to the 
gold surface, our fc-point sampling is equivalent to having 
(6-2)x(6-2) = 12xl2 fc-points in a lxlxl real-space gold 
cell. A further increase in the number of fc-points does 
not affect our results. 

In a self-consistent calculation on the whole structure, 
the first-order Methfessel-Paxton scheme is used with 
the smearing parameter er=6.0 meV. This choice of a 
is necessary to resolve individual energy levels in pro- 
jected density of states (PDOS) and to compute the MAB 
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caused by the SOI in the Mn.12 molecule. The total mag- 
netic moment is set to 20 [is and the discussion in Sec. Ill 
suggests that 20 [Ib is the ground-state total magnetic 
moment for the whole structure. The combination of 
high energy cutoffs and a large number of valence elec- 
trons involved indicates the numerically expensive nature 
of the performed calculations so we do not further relax 
the whole structure. A comprehensive analysis of the 
forces exerted on the structure is provided. The average 
over absolute values of all of the x, y, and z force compo- 
nents is (0.03,0.03,0.05) eV/A with standard deviation 
of (0.08,0.07,0.08) eV/A. Out of a total of 304 atoms, 
forces on only 25 atoms are greater than 0.15 eV/A and 
the large forces occur near the sulfur atoms. The x and y 
components of the forces on the gold atoms bonded to the 
sulfur atoms are 0.60 eV/A and the z force components 
are 0.20 eV/A with alternating signs. The forces on near- 
est neighbors to those gold atoms are 0.20 eV/A along 
the x and y axes. The force on one of the sulfur atoms has 
components (0.57,0.29,0.13) eV/A because of the mis- 
match of placing the sulfur atom at a hollow site of the 
gold slab. The forces onto the two carbon atoms bonded 
to the sulfur atoms are 0.39 eV/A along the z axis. The 
oxygen atoms bonded to those carbon atoms have z force 
components that fall between —0.66 and —0.52 eV/A. 
The Mn ions connected to the S via C and O atoms (two 
of the Mn(I) and two of the Mn(II) ions) have z force 
components of 0.24 eV/A. The analysis of the forces 
indicates that relaxation of the whole structure may pos- 
sibly break the symmetries of the Mn(I) and Mn(II) ions, 
and could further modify its magnetic properties. 



III. RESULTS AND DISCUSSION 

The electronic structures of the six gold monolayers, 
geometries 1 and 2, and the whole structure are cal- 
culated and discussed using the PDOS onto particular 
orbitals of specific atoms. A comparison of the PDOS 
for the whole structure (Fig. 0]) with those for geometry 
2 (Fig. 0, computed under identical conditions, reveals 
that the molecule of geometry 2 is weakly coupled to the 
gold surface for the given link, and that the molecular or- 
bitals broaden by at most an order of 0.1 eV. The Fermi 
level of the gold slab is well above the majority LUMO 
but below the minority LUMO of geometry 2 (Fig. [3]) . As 
a result, only electrons with majority spin can be trans- 
ferred from the gold slab to the molecule. The ground- 
state total magnetic moment can also change upon ad- 
sorption on the surface due to strong bonds formed be- 
tween the sulfur atoms and the gold slab. The amount 
of the charge transfer is quantified from two methods: 
(i) a onc-dimcnsional electronic charge density difference 
between the whole structure and geometry 2, and (ii) a 
shift in the lowest-lying molecular orbital as a function of 
extra partial charges. The change in the total spin mag- 
netic moment is obtained from a one-dimensional spin 
density difference. The effect of the charge transfer and 
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FIG. 4: (Color online) Projected densities of states onto the 
majority and minority Mn d, O p, and S p orbitals for the 
whole structure: geometry 2 on Au(lll). The zero in the 
horizontal scale denotes the Fermi level. 
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FIG. 5: (Color online) Projected densities of states onto the 
majority and minority Mn d, O p, and S p orbitals for ge- 
ometry 2. The zero in the horizontal scale denotes the mid- 
point between the HOMO and LUMO energies. The ovals and 
squares correspond to the Mn d and O p orbitals hybridized 
with the S p orbitals. Notice that the vertical scale for the S 
orbitals here differs from that in Fig(4] 



spin magnetic moment change on the MAB is discussed 
from self-consistent SOI calculations. 



A. Projected densities of states (PDOS) 

The PDOS shown in Figs. QJ [5] and [S] are generated with 
100 data points/eV and smearing parameter a=6 meV 
without the SOI, over atom-centered spheres with the 
radii given in Table |TJ The horizontal scale in Figs. [5] 
and [5] is shifted by 0.34 eV and 3.5 meV, respectively, 
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FIG. 6: (Color online) Projected densities of states onto the 
majority and minority Mn d and O p orbitals for geometry 1. 
The zero in the horizontal scale denotes the midpoint between 
the HOMO and LUMO energies. 



such that the zero in the horizontal scale corresponds 
to the midpoint of the HOMO and LUMO energies. The 
PDOS onto the gold s and d orbitals (not shown) indicate 
a large density of states within 0.5 eV below and above 
the gold Fermi level. To directly compare the PDOS for 
the whole structure with those for geometries 1 and 2, 
exactly the same parameter values are used for the three 
systems, such as an identical unit cell, fc-point sampling, 
energy cutoff, cutoff in the augmentation charges, PDOS 
sampling parameters, smearing parameter, and total en- 
ergy convergence tolerance. For geometry 2 the HOMO 
and orbitals right below the HOMO are mostly from the 
Mn(III) ions, neighboring O anions, and S atoms, while 
the LUMO and orbitals right above the LUMO are from 
the Mn(II) ions, neighboring O anions, and the S atoms. 
As highlighted by ovals and squares in Fig. the S p 
orbitals are hybridized with the Mn d and O p orbitals. 
Within 0.5 eV below the midpoint of the HOMO and 
LUMO energies, the spin moments of the S atoms are 
antiparallel to those for the Mn(III) and O atoms. For ge- 
ometry 1 the HOMO and orbitals right below the HOMO 
are mainly from the Mn(II) ions and O anions, while the 
LUMO and orbitals right above the LUMO are from the 
Mn(I) and Mn(III) ions and O anions (Fig. [BJ. These 
main differences of the PDOS between geometries 1 and 
2 are linked to the difference of the total ground-state 
magnetic moment between them. 

For the whole structure the PDOS (Fig. 0| onto the 
Mn d and O p orbitals noticeably differ from those for 
geometry 2 because of the difference in the total mag- 
netic moment between them. Although the PDOS for the 
whole structure are similar to those for geometry 1, the 
peak heights in the former are reduced and the orbitals 
broaden compared to the latter. The broadening is of an 
order of 0.1 eV (Figs. [31 [U and [5]). The single-electron 
charging energy for geometry 2 is calculated from the 
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FIG. 7: (Color online) Red filled circles: Charge density dif- 
ference between the whole structure and geometry 2 vs the z 
coordinate. A, B, and C denote the locations of the S atoms 
(at 2=0.09 A), the bonding C atoms, and the Mn ions illus- 
trated in the inset of Fig. [9] Upon integration from z — 0.09 
(z = -0.93) to 2=12.00 A, a total of 1.23 (23.15) electrons 
is transferred from the surface to the molecule. Black empty 
squares: Charge density difference between geometries 3 and 
2. Blue dashed curve: Charge density difference between a 
molecule with extra 0.3 electrons added to geometry 2 (de- 
noted geometry 4) and a neutral molecule with geometry 2. 



energy difference between a neutral molecule and one 
with an extra-electron added. When one extra electron 
is added to geometry 2, the total ground-state magnetic 
moment increases to 19 (is- Considering the change in 
the total magnetic moment, we obtain the charging en- 
ergy of 3.7 eV. This energy is one order of magnitude 
greater than the estimated broadening of the molecular 
orbitals so the molecule is weakly bonded to the surface. 



B. Charge transfer 

The charge transfer between the surface and the 
molecule of geometry 2 is computed from a one- 
dimensional electronic charge density difference between 
the whole structure and geometry 2 along the z axis. 
The positive sign in the charge difference means an elec- 
tronic charge is transferred from the gold surface to the 
molecule. The charge difference shown as the curve with 
filled circles in FigJTl increases abruptly towards the sul- 
fur atoms (vertical line A) due to the electronic contri- 
bution of the 36 gold surface atoms as well as the sul- 
fur atoms. (Each gold atoms has 11 valence electrons.) 
The deep penetration of the charge of the gold atoms 
beyond the sulfur atoms makes it difficult to precisely 
calculate the amount of the charge transfer. Including 
the long tail of the charge of the gold surface atoms, a 
charge of 1.23 electrons (23.15 electrons) is transferred 
from the surface to the molecule when the integration 
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An (e) 

FIG. 8: (Color online) Shift of the lowest-lying molecular or- 
bital (MOl) energy for geometry 2, AMOl, as a function of 
extra partial electronic charge, An. The thick horizontal line 
indicates the energy difference of the MOl between the whole 
structure and geometry 2. 



is performed from the sulfur atoms, z — 0.09 (the mid- 
distance between the first gold layer and sulfur atoms, 
2 = -0.93) to 2=12.00 A (Fig. EJ). It is found that the 
transferred charge is mostly localized in the sulfur and 
bonding carbon and oxygen atoms, while little charge is 
transferred to the Mn ions. A similar trend is shown for 
the charge density difference between geometry 2 and the 
structure where two H atoms are bonded to the S atoms 
[denoted geometry 3, Fig.OJb)]. In this case, the charge 
difference obtained from integration from z = 0.09 to 
z = 12.00 A, is 0.17 electrons. Notice that these two 
sets of the charge density difference match well above 
vertical line B in Fig. \7\ indicating that the electronic 
environment around the molecular magnetic core is very 
similar in the whole structure and geometry 3. To further 
stress this point, we calculate the charge density differ- 
ence between a neutral molecule with geometry 2 and the 
same molecule with 0.3 free electrons added. The extra 
free electrons are uniformly distributed over all atoms, as 
shown in Fig. 

The charge difference is also calculated from the change 
of the lowest- lying occupied molecular orbital (MOl) of 
geometry 2 as a function of extra partial electrons added 
to the molecule, following the procedure described in 
Ref. H3. As shown in Fig. [51 the MOl changes linearly 
with extra partial charges. Here a correction due to the 
discontinuity in the derivative of the exchange-correlation 
potential with respect to electron density within DFT is 
not taken into account. The molecular orbitals of geome- 
try 2 are modified upon adsorption onto the gold surface 
due to the interaction with the surface. Lying deeper in 
energy than the first occupied orbital of the gold slab, 
the MOl is not hybridized with the gold surface. Thus, 
the shift of the MOl in energy upon adsorption on the 
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FIG. 9: (Color online) Red filled circles: Spin magnetic mo- 
ment difference between the whole structure and geometry 2. 
The difference of the total magnetic moment is 1.80pB when 
the spin density difference is integrated from z = —0.93 to 
2=12 A. Black empty squares: Magnetic moment difference 
between geometry 3 and 2. Blue dashed curve: Magnetic mo- 
ment difference between geometries 4 and 2. The inset shows 
lines A and B and region C. 



surface can be considered purely due to the charge trans- 
fer from the surface to the molecule, assuming that the 
Coulomb repulsion is small. The calculated shift of the 
MOl between the whole structure and neutral geometry 
2, 363 meV, leads to a charge transfer of 0.28 electrons, 
as depicted in Fig. [SJ 



C. Spin magnetic moment difference 

The spin magnetic moment difference between the 
whole structure and geometry 2 is shown in Fig. O The 
increase in the total magnetic moment of the molecule 
upon adsorption amounts to 1.80 fis when the one- 
dimensional spin density difference is integrated from 
z = -0.93 to 2=12 A. As shown in Fig. M and Table HH 
the major difference arises from the sulfur atoms due to 
the strong bonds to the gold surface, while the next con- 
siderable difference occurs at the O anions and Mn ions 
close to the surface. Since only electrons with majority 
spin can be transferred to the molecule, the Mn(I) and 
Mn(II) ions close to the surface have higher magnetic 
moments than those for geometry 2. The calculated dif- 
ference of the spin magnetic moment corroborates the 
difference of the total ground-state magnetic moment be- 
tween the whole structure and geometry 2. It is empha- 
sized that this large change of the total magnetic moment 
is mainly caused by rearrangement of the magnetic mo- 
ments of the S atoms and Mn ions due to the strong bonds 
to the gold surface, rather than by the charge transfer. 
Note that the total electronic charge must be conserved 
but not necessarily the total spin magnetic moment. A 
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TABLE III: Magnetic anisotropy barrier obtained from self- 
consistent SOI calculations for isolated Mni2 molecules (ge- 
ometries 1, 2 and 3) and the whole structure, in units of K. 



Geometry 1 


Geometry 2 


Geometry 3, whole structure 


66.7 


66.9 


60.7 



tal magnetic moment of 20 hb, the calculated MAB for 
geometry 3 is 60.7 K that is 6 K lower than the MAB for 
geometries 1 and 2 (Table [TTT|) . Thus, a reasonable value 
for the whole structure would be also 60.7 K. The reduc- 
tion in the barrier is due to the decrease in the single-ion 
anisotropy for the two Mn(II) ions close to the surface. 
IV. CONCLUSION 



spin magnetic moment of 0.20 [is is obtained for the 
gold slab up to z = —0.93 A. As illustrated in Fig. [5J 
the magnetic moment difference between geometries 3 
and 2, is remarkably similar to that between the whole 
structure and geometry 2. The two sets of the magnetic 
moment differences are, however, qualitatively dissimi- 
lar to the magnetic moment difference caused by 0.3 free 
electrons, shown by the dashed curve in Fig. [5J The ex- 
treme similarity in the one-dimensional charge and mag- 
netic moment distributions in geometry 3 and the whole 
structure, as well as the moderate level broadening in the 
presence of gold, suggest that one estimate the MAB for 
the whole structure from a calculation on geometry 3. 

D. Magnetic anisotropy barrier (MAB) 

The magnetic anisotropy comes from Jahn- Teller dis- 
tortion around the eight Mn 3+ ions at the outer ring. It 
is known experimentally that for geometry 1 the second- 
order magnetic anisotropy contributes to the total MAB 
by 55 K and the fourth-order anisotropy by 10 KJ£ 
With the total magnetic moment fixed, the MAB is com- 
puted from the change of the total energy upon a ro- 
tation of the spin quantization axis from the z to the 
x axis. The MAB calculated from a non self-consistent 
SOI calculation is 57.2 K for geometry 1, while the bar- 
rier from a self-consistent SOI calculation increases up to 
66.7 K. The latter is more precise and may include higher- 
order magnetic anisotropy. This result agrees with the 
all-electron DFT-calculated barrier— and experimental 
values. 18,36 Henceforth we report only the MAB obtained 
self-consistently. The MAB with total magnetic moment 
of 18 [is is 66.9 K for geometry 2 so it remains the 
same as that for geometry 1. The transverse magnetic 
anisotropy is calculated to be 0.02 K, which is caused by 
the sulfur atoms. As discussed, geometry 3 is used to 
calculate the MAB for the whole structure. With the to- 



We have investigated, within DFT, the interaction be- 
tween a single-molecule magnet Mni2 monolayer and a 
Au(lll) surface through thiol bonds without long alkane 
chains. Despite a very short bond length between the 
sulfur-terminated Mni2 molecule and the gold slab, the 
broadening of the molecular orbitals was much less than 
the single-electron charging energy of the molecule. In 
the ground state the sulfur-terminated Mni2 molecule 
has a total magnetic moment of 18 /is, but its total mag- 
netic moment increases to 20/is upon adsorption onto 
the gold surface. The noticeable charge transfer from 
the surface to the Mni2 is attributed to (i) the relative 
position of the gold Fermi level to the LUMO and (ii) 
the long tail in the electronic cloud at the boundary of 
the gold surface. The self-consistent SOI calculation sug- 
gested a decrease in the MAB by 9% and a significant 
transverse magnetic anisotropy for the whole structure 
compared to the isolated Mni2. Although a relaxation of 
the whole structure may bring additional transverse mag- 
netic anisotropy caused by minor symmetry breaking of 
the Mn ions close to the surface, it is unlikely that the 
molecules are completely collapsed upon adsorption. The 
analysis and results given in this study are not limited 
to the Mni2 molecules and may be applicable to systems 
where other types of magnetic molecules or SMMs are 
adsorbed onto nonmagnetic surfaces or bridged between 
nonmagnetic electrodes. 
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